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We show that in suitable anisotropic ferromagnets, both stable and metastable "tilted" phases 
occur, in which the magnetization M makes an angle between zero and 180 degrees with the exter- 
nally applied H. Tuning either the magnitude of the external field or the temperature can lead to 
continuous transitions between these states. A unique feature is that one of these transitions is be- 
tween two metastable states. Near the transitions the longitudinal susceptibility becomes anomalous 
with an exponent which has an exact scaling relation with the critical exponents. 
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Anisotropic anti- ferromagnets 0, have long been 
known to have very rich phase diagrams, including phases 
exhibiting Ising, XY, and both types of order simultane- 
ously. 

Less attention seems to have been focussed on the fer- 
romagnetic case. Here we analyze that case, and show 
that it, too, can exhibit phases with simultaneous Ising 
and XY order. In this case, such "mixed" order implies 
extremely novel "tilted" phases, in which the magneti- 
zation M makes an angle between zero and 180 degrees 
with the externally applied H. These phases can be both 
stable and metastable. 

Our results are summarized in Fig. [T] which illustrates 
both the equilibrium and metastable phases of a suitable 
hexagonal ferromagnetic crystal in the H-t± plane, where 
H is the external magnetic field, and ij_ is a phenomeno- 
logical parameter that increases monotonically with tem- 
perature T. "Suitable" , in this context, means the crystal 
field obeys certain conditions that we will specify more 
precisely later. For now, we simply point out that these 
conditions prove to be "generic" : that is, they do not 
require "fine-tuning" of any material parameters. This 
does not mean that all hexagonal crystals will exhibit 
the phase diagram Fig. [TJ it simply means that some of 
them should. 

Metastable phases, of course, depend on the sample 
history. In Fig. [TJ we have assumed that the system 
starts in an ordered state with both the external field 
H and the magnetization M large, and in the —z direc- 
tion, where z is the unit vector normal to the hexagonal 
planes. Keeping temperature fixed, and the external field 
H along the z axis, the component H z of H along z is 
then varied from large negative to large positive values. 

In the region below the locus FDG and above the locus 
EBKAJ, there is no metastability, and the equilibrium 
state is just the conventional one, with the magnetization 
M || H, the external field. The region JAKDG is like- 
wise quite familiar: here, the equilibrium state remains 
M II H, while the metastable state is the conventional 



one with M || -H. 

All of the other regions of the phase diagram exhibit 
tilted phases either in the equilibrium or the metastable 
state, or both. In the region BKD, the tilted phase is 
metastable; the equilibrium state still has M j H. In 
the regions CDF and CBE, there is no metastability, 
the equilibrium states are tilted; while in CBD, both 
the metastable and equilibrium states are tilted, but the 
components of M along the external field are of opposite 
sign! 

The loeii DF and EB are equilibrium tilting transi- 
tions, while the locus KD is a very strange one indeed: 
it represents a purely metastable transition: i.e., a con- 
tinuous transition between two metastable states. 

For the sample history we specified, the locus BD 
is not a transition line if the system is trapped in the 
metastable states. However, if the system is allowed to 
equilibrate, it is also a equilibrium transition. 

In all of the tilted states, the projection of the tilted 
magnetization onto the hexagonal planes always lies 
along one of six six- fold symmetry related directions. 

For a part of the parameter space of our model, the 
tilting is continuous, and belongs to the universality class 
of the three-dimensional XY model[|[; that is, for H —> 
H c (t±), where H c (t±) is the i^-dependent critical field H 
at which the tilting transition occurs (i.e., value of H on 
one of the tilting loeii DF, EB, and KD just discussed), 
the tilt angle 9 is given by 

e~\H-H c f, (i) 

where the universal exponent j3 = 0.3485 ± 0.00020 is 
the order parameter exponent for the three-dimensional 
XY model. 

The tilting transition can also be crossed by varying 
temperature. In this case, H and H c are replaced in 
Eq. fTJ) with T and T c , respectively, with T c being the 
temperature on the tilting loch. Note that, in contrast 
to the usual ferromagnet-paramagnet critical point, tem- 
perature T and external field H are equivalent here, in 
the sense just described. 
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Since the projection M± of the magnetization M per- 
pendicular to the applied field is the order parameter for 
this tilting transition, both the associated susceptibility 
(namely, the uniform transverse susceptibility) \±- an d 
the correlation length £ for correlations of M± diverge 
near the tilting transition, according to the laws: 



f ~ \H - H c 



(2) 



where the critical exponents 7 = 1.3177 ± 0.0005 and 
v = .67155 ± 0.00027 are respectively the universal sus- 
ceptibility and correlation length exponents of the three- 
dimensional XY modelQ. 

In addition, the generalized longitudinal susceptibility 
is renormalized by the critical fluctuations and becomes 
wavelength-dependent. It displays a weak anomaly 



Xz(q) 



constant 



C ± \H-H C 



9»r 1 f31 



where a = —0.0146 ±0.0008 is the universal specific heat 
exponent of the three-dimensional XY model. C± are 
non-universal positive constants for H above and below 
H c , respectively. Their ratio jt- is universal. 

We have also investigated these phenomena for other 
crystal lattice symmetries. For a cubic ferromagnetic 
crystal we find that if the external field is along one of 
the cubic axes (e.g., (100)), a continuous tilting transi- 
tion is possible. The universality class of the transition is 
again that of the three-dimensional XY model, and the 
longitudinal susceptibility is again given by ©. 

If the external field is along one of the body diagonal 
directions (e.g., (Ill)), the tilting transition is in the uni- 
versality class of the three-dimensional three-state Potts 
model, which is believed Q to be first order. 

In the case of an orthorhombic crystal, when the exter- 
nal field is along one of the three non-equal primary axes, 
the tilting transition can be continuous. Furthermore, it 
can happen between metastable as well as equilibrium 
states. The universality class of this transition is that 
of the three-dimensional Ising model. Near the critical 
point Xz displays a divergence of the form , but with 
C± < 0, v = 0.630 ± 0.001, and a = 0.109 ± 0.004@, 
where v and a are respectively the universal correla- 
tion length and specific heat exponents of the three- 
dimensional Ising model. 

Our model for an anisotropic hexagonal ferromagnet 
in an external field H is: 



F 



1 



d d r 



g z M 2 + t±\M ± \ 2 + C\VM\ 2 + u z M 4 z 



+u±\M±\ 4 + u± z \M±\ 2 M 2 - 2HM Z 



(4) 



where _L denotes the x — y plane, taken to be the plane 
of hexagonal symmetry, and the positive direction of the 
field is along z. This model with H = was first studied 
by Fisher et.al.[l| as a model for antiferromagnets in uni- 
axial crystals. It includes all terms to fourth order in M 




FIG. 1: Phase diagram for hexagonal ferromagnets below the 
Curie temperature in the presence of an external field H per- 
pendicular to the hexagonal planes (this direction will here- 
after be called the "z-axis"). Each phase is identified by two 
arrows. The right arrow denotes the orientation of the magne- 
tization in the true equilibrium state, and the left one denotes 
the orientation in the state actually reached upon increasing 
H z from large negative values; when this differs from the true 
equilibrium state, the state is metastable. The locus KD is 
a "meta-critical" line separating two metastable states, while 
the loeii BE and DF are "true" critical lines separating dis- 
tinct equilibrium states. 



allowed by the hexagonal symmetry. There are sixth or- 
der in M terms allowed by the hexagonal symmetry that 
break the continuous rotation invariancc of this model in 
the _L-plane down to six-fold rotational invariance. Such 
terms pick out six equivalent preferred directions within 
the _L-plane for the tilting, but do not affect either the 
topology of the phase diagram in the figure, or the uni- 
versality classes of any of the various transitions therein. 

In what follows we will consider only the case g z < 
0. Moreover, for simplicity we restrict ourselves to the 
region of parameter space u± tZ > 0, < u± z < 2^/u±u z . 

We take the initial field to be so strong (i.e., a large 
negative H) that in the ground state the nonzero mag- 
netization points along —z. 

We will begin by treating this model in Landau theory, 
in which we find the state of the system by minimizing 
this Landau free energy Eq. Expanding it around 

the ground state by writing M = (—M + 5M z )z + Mj_ 
where M satisfies M g z + 2u z M$ + H = 0, we obtain 



A{5M Z ) 2 + B(SM Z )\M ± \ 2 + C7|Vi\/ ± | 2 
+D\M ± \ 2 +u ± \M ± \ 4 } (5) 



where we have defined A = g z + 6MqU z , B = 2Mqu± z , 
and D = t± + u± z M 2 . Initially, (i.e., when H is large 
and negative), we have A, D > 0. 

Now let us consider increasing the value of H gradually 
while fixing all other parameters (e.g., the temperature). 
Increasing H decreases Mo, hence decreasing the coeffi- 
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cient D. We find that for t± < 0, when H reaches 



-(Jr. 



2uz 

U±2 



-t± 



-t± 



(6) 



D = 0. H = H c defines the tilting locus AKDF in the 
figure. The locus DE is the mirror image of DF about 
the t± axis. For H > H c , D < 0, and the magnetization 
tilts away from — z such that its transverse component 
Mj_ becomes nonzero. For small tilting, the tilt angle 
6 oc \M±\/Mq. Therefore, M± is the order parameter of 
this transition, and 9 is proportional to its magnitude. 

However, if the coefficient A changes sign before D 
does, the untilted phase becomes unstable. This hap- 
pens at H = Hf = -a^J^al, Thus for the tilting 

transition to occur we must have H c < Hf. In terms of 
the parameters in model ([Sjl. this condition can be writ- 
ten as t± < t max = ^fg z - For t± > t max , increasing H 
only leads to the complete flipping of the magnetization 
from — z to z without any tilting. This flipping occurs on 
the locus A J in the figure. Assuming t± < t max , near the 
tilting transition we can eliminate the degree of freedom 
5M Z from the free energy Eq. ([5]). This gives: 



F 



d d r 



C\VM ± \ 2 +D\M A 



+ u± 



AA 



IMi 



(7) 



where D ~ H c — H near the transition. For the tilting 
transition to be continuous, at H = H c the coefficient of 
the quartic term has to be positive, which leads to 



t± < 



g z u±u± z 



(8) 



This condition can be satisfied in the region on the left 
side of the point K in the figure. Thus we conclude that 
a continuous tilting transition is possible on the locus 
KDF. Since M± is a two-component vector, and since 
the model Eq. Q is invariant under rotations of M± in 
its plane, the universality class of the tilting transi- 
tion is that of the three-dimensional XY model. Recog- 
nizing this leads to the critical exponents for this transi- 
tion quoted in Eqs. (JTJ, ©, and © at the beginning of 
this paper. 

Further increasing H after the continuous tilting tran- 
sition eventually leads to an instability of the tilted 
phase. The determinant of the Hessian matrix of the 
Landau free energy Eq. (j4]) vanishes at this instability 
limit, which, after some algebra, we find is at 



H 



1 



U±y^Au±U z — U 



2 

±z 



t±uj_ z - 2g z u± 



(9) 



This instability limit is illustrated by the locus KBC in 
the figure. 

Our main results are summarized by the phase dia- 
gram in H-t±_ plane illustrated in the figure. The other 



features of this phase diagram, including the locii of all 
the phase boundaries and special points, can be obtained 
by straightforward minimization of the Landau free en- 
ergy. Details will be given in a future publication Q. 

Since the lifetimes of the metastable states are finite, 
the transitions between them will be observable only if 
the experiment is conducted within their lifetimes. 

Because the Landau coefficients g z , t±, u z , u± z are 
all temperature dependent, the tilting transition can be 
induced by tuning the temperature as well. For conve- 
nience, in what follows we consider the tilting transition 
induced by fixing H and varying T only. The correspond- 
ing results for the opposite case (i.e., fixing T and varying 
H) can be derived in essentially the same way. 

Fluctuations of 5M Z arc greatly affected near the tilt- 
ing transition by the critical fluctuations of M±. To see 
this, note that the model Eq. ([S]) implies that the truly 
massive quantity is the combination SM Z + B\M±\ 2 /2A. 
Therefore, 



(SM z ) = -^( 



constant 



(10) 



where the ± subscript distinguishes coefficients above 
and below T c , respectively. These coefficients are dif- 
ferent, and non-universal; however, their ratio is 
universal, and identical to the analogous ratio of the spe- 
cific heat coefficients above and below T c . 

In deriving the above equation we have used the well- 
known result 0, from the theory of critical phenomena: 



A T± \T-T C 



constant . 



(11) 



and defined A L± = -BA T± /2A. 

Now we can calculate the average value of M z using 
(M z ) = -M + (8M Z ). This gives 



(M z 



-M (H,T) + A L ±\T-T C 



l-a 



constant. (12) 



Thus we see that M z exhibits a singularity at the tilt- 
ing transition, despite the fact that M z itself is not the 
order parameter for this transition. Indeed, this sin- 
gularity (|T2"j) of M z is very similar to, and arises from 
the same mechanism as, the well- known [|| singularity of 
the lattice constant in a compressible ferromagnet at the 
ferromagnet-paramagnet transition^ . 

Right at the critical point, where D = 0, the model 
Eq. ((Sjl becomes very similar to that for an isotropic 
ferromagnet ordered with M \\ z. In such an isotropic 
ferromagnet, the longitudinal susceptibility \z is diver- 
gently renormalized by the transverse fluctuations [lfl] . 
We expect a similar anomaly here, with the crucial dif- 
ference that, in the isotropic case, rotational invariance 
requires that the coefficients in model Eq. (|5|) satisfy 
A = B = Au±_ , while in our problem, B and u± are inde- 
pendent, since rotational invariance is broken even when 
D — 0. Similar issues arise in the smectic-A-smectic-C 
phase transition in anisotropic environments [ill [l2j . 
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We calculate this anomaly as follows: In real space 
Xz{r,r') gives the change in (M z (r)} in response to an 
external field which is along z and at r ' . That is, 

TXz(r,f') = ([SM z (f) - (SM z (r))][5M z (f') - 

(SM z (r'))]). (13) 

Noting that the combination 5M Z + is massive, 

we write SM. 



B\M 



Xz(q) 



2^4 — ; then, in Fourier space 
D 2 



AA 2 T 



G{q), 



(14) 



where G(q) is the Fourier transform of the correlation 
function {jM ± (r)| 2 |ML(0)| 2 ^ - ^|A/_l| 2 ^ . Standard 
scaling arguments Q imply that near the transition 



(15) 



where n = — a/V, and C'± = At±. Plugging (fi"5|) into 
(HU) leads to the result ©. 



The results described earlier for the tilting transition 
in cubic and orthorhombic ferromagnets can be derived 
by applying the type of analysis used above for hexagonal 
ferromagnets to Landau theories that respect the (differ- 
ent) symmetries for those cases. This will be discussed 
in a future publication Q. 

In summary, we have studied the field induced phase 
transitions between ordering states in anisotropic hexag- 
onal ferromagnets. Some of these transitions are sig- 
naled by a tilting of the orientation of the magnetiza- 
tion; in some cases these tilting transitions occur between 
metastable states. We also found the universality classes 
of these transitions. Similar scenarios can occur for a 
variety of crystal symmetries and field orientations Q- 

The ideas presented here have implications for fer- 
romagnetic superconductors [13, EH which exhibit a 
spontaneous-magnctization-induccd flux lattice. Since 
these materials are anisotropic ferromagnets, they should 
also undergo a tilting transition in an applied external 
field. The coupling between the magnetization and the 
flux lines will change the tilting transition. We leave this 
interesting problem for future workQ- 



[1] J. M. Kosterlitz, David R. Nelson, and Michael E. Fisher, 

Phys. Rev. B. 13, 412(1976). 
[2] A. D. Bruce and A. Aharony, Phys. Rev. B. 11 478 (1975) 
[3] For a discussion of universality classes, see, e.g., P. M. 

Chaikin and T. C. Lubensky, Principles of Condensed 

Matter Physics, Cambridge (1997). 
[4] M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, 

and E. Vicari, Phys. Rev. B 63, 214503 (2001). 
[5] D. Mukamel, M. E. Fisher, and E. Domany, Phys. Rev. 

Lett. 37, 565 (1976). 
[6] R. Guida and J. Zinn- Justin, J. Phys. A 31, 8103 (1998). 
[7] L. Chen and J. Toner, in preparation. 
[8] See, e.g., A. T. Dorsey, P. M. Goldbart, and J. Toner, 

Phys. Rev. Lett. 96, 055301 (2006). 
[9] D. J. Bergman and B. I. Halperin, Phys. Rev. B 13, 2145 

(1975). 

[10] G. F. Mazenko, Phys. Rev. B 14, 3933 (1976). 

[11] G. Grinstein and R. A. Pelcovits, Phys. Rev. A 26, 2196 



(1982). 

[12] L. Chen and J. Toner, Phys. Rev. Lett. 94, 137803 
(2005); L. Chen and J. Toner, Phys. Rev. E. 79, 031703 
(2009). 

[13] A. Huxley, I. Sheikin, E. Ressouche, N. Kernavanois, D. 
Braithwaite, R. Calemczuk, and J. Flouquet, Phys. Rev. 
B 63, 144519 (2001); S. S. Saxenaet al., Nature (London) 
406, 587 (2000); D. Aoki, A. Huxley, E. Ressouche, D. 
Braithwaite, J. Floquet, J. P. Brison, E. Lhotel, and C. 
Paulsen, Nature (London) 413, 613 (2001); C. Pfleiderer, 
M. Uhlarz, S.M. Hayden, R. Vollmer, H. von Lohneysen, 
N. R. Bernhoeft, and G.G. Lonzarich, Nature (London) 
412, 58 (2001) 

[14] L. Radzihovsky, A. M. Ettouhami, K. Saunders, and 
J. Toner, Phys. Rev. Lett. 87, 27001 (2001); A. M. 
Ettouhami, K. Saunders, L. Radzihovsky,and J. Toner, 
Phys. Rev. B 71, 224506 (2005). 



